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Abstract 

The Lagrangian-Hamiltonian unified formalism of R. Skinner and R. Rusk was originally 
stated for autonomous dynamical systems in classical mechanics. It has been generalized for 
non-autonomous first-order mechanical systems, as well as for first-order and higher-order 
field theories. However, a complete generalization to higher-order mechanical systems has 
yet to be described. In this work, after reviewing the natural geometrical setting and the 
Lagrangian and Hamiltonian formalisms for higher-order autonomous mechanical systems, 
we develop a complete generalization of the Lagrangian-Hamiltonian unified formalism for 
these kinds of systems, and we use it to analyze some physical models from this new point 
of view. 
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1 Introduction 

In recent decades, a strong development in the intrinsic study of a wide variety of topics in 
theoretical physics, control theory and applied mathematics has been done, using methods of 
differential geometry. Thus, the intrinsic formulation of Lagrangian and Hamilonian formalisms 
has been developed both for autonomous and non-autonomous systems. This study has been 
carried out mainly for first-order dynamical systems; that is, those whose Lagrangian or Hamil- 
tonian functions depend on the generalized coordinates of position and velocity (or momentum). 
From the geometric point of view, this means that the phase space of the system is in most cases 
the tangent or cotangent bundle of the smooth manifold representing the configuration space. 
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However, there are a significant number of relevant systems in whicli the dynamics have ex- 
plicit dependence on accelerations or higher-order derivatives of the generalized coordinates of 
position. These systems, usually called higher-order dynamical systems, can be modeled geomet- 
rically using higher-order tangent bundles [17j. These models are typical of theoretical physics; 
for example those describing the interaction of relativistic particles with spin, string theories 
from Polyakov and others, Hilbert's Lagrangian for gravitation or Podolsky's generalization of 
electromagnetism (see [8] and references cited there). They also appear in a natural way in 
numerical models arising from the discretization of first-order dynamical systems that preserve 
their inherent geometric structures |16j . There are a lot of works devoted to the development of 
the formalism of these kinds of theories and their application to many models in mechanics and 
field theory (see, for instance, [2], [3], [1], [7], [9], [E], 120], [25], [26], [35], [37], [38]). 

Furthermore, a generalization of the Lagrangian and Hamiltonian formalisms exists that 
compresses them into a single formalism. This is the so-called Lagrangian-Hamiltonian unified 
formalism, or Skinner-Rusk formalism due to the authors' names of the original paper. It was 
originally developed for first-order autonomous mechanical systems [39], and later generalized to 
non-autonomous dynamical systems [6l , control systems [5] , first-order classical field theories 
[151118] and, more recently, to higher-order classical field theories [101141]. Nevertheless, although 
the geometrization of both higher-order Lagrangian and Hamiltonian formalisms was already 
developed for autonomous mechanical systems [111 |T71 [23l |29] , a complete generalization of the 
Skinner-Rusk formalism for higher-order mechanical systems has yet to be developed. A first 
attempt was outlined in |13] . with the aim of providing a geometric model for studying optimal 
control of underactuated systems, although a deep analysis of the model and its relation with 
the standard Lagrangian and Hamiltonian formalisms was not performed. 

Thus, the aim of this work is to provide a detailed and complete description of the Lagrangian- 
Hamiltonian unified formalism for higher-order autonomous mechanical systems. Our approach 
is different from that given in p!3] (these differences are commented on Section [5]) . 

The paper is organized as follows: 

Section [2] consists of a review of the basic definitions and the geometric structures of higher- 
order tangent bundles, some of which are generalizations of the geometric structures of tangent 
bundles; namely, the canonical vector fields, the almost-tangent structures and semisprays; 
whereas others such as the Tulczyjew derivation are needed for developing the Lagrangian and 
Hamiltonian formalisms of higher-order mechanical systems, which are also described in this 
section. In particular, higher-order regular and singular systems are distinguished. 

The main contribution of the work is found in Section [3l where the geometric formulation of 
the Lagrangian-Hamiltonian unified formalism for higher-order autonomous mechanical systems 
is described in detail, including the study of how the Lagrangian and Hamiltonian formalisms 
are recovered from that formalism. 

Finally, in Section HI two examples are analyzed in order to show the application of the 
formalism; the first is a regular system, the so-called Pais-Uhlenbeck oscillator, while the second 
is a singular one, the second-order relativistic particle. 

The paper concludes in Section [5] with a summary of results, discussion and future research. 

All the manifolds, the maps and the structures are smooth. In addition, all the dynamical 
systems considered are autonomous. Summation over crossed repeated indices is understood, 
although on some occasions the symbol of summation is written explicitly in order to avoid 
confusion. 
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2 Higher-order dynamical systems 

2.1 Geometric structures of higher-order tangent bundles 
(See nil [Ml [23 [Ml [2Z] for details). 

2.1.1 Higher-order tangent bundles 

Let Q be a n-dimensional differentiable manifold, and A; G N. The kth-order tangent bundle 
of Q, denoted by T^Q, is the {k + 1) n-dimensional manifold made of the /c-jets with source at 
G M and target Q; that is, T'^Q = J^{R, Q). It is a submanifold of J^{R, Q). 

We have the following canonical projections: if r < A;, 

T^Q T'^g , /?'=: T^Q Q 

where (5"'^(0) denotes a point in T^Q; that is, the equivalence class of a curve cr : / C M — )• Q by 
the A: -jet equivalence relation. Notice that Pq = j3^, where T^Q is canonically identified with Q, 
and = Id^kQ. 

If {U,ip) is a local chart in Q, with ip = {y^^), 1 < ^4 < n, and o": M — ?• Q is a curve in Q 
such that (t(0) G U; by writing cr^ = 93^ o fj, the A-jet d-^{0) is given in {f3^)~^ (U) = T^U by 

(jig- A 

((7^, (7^, . . . , g^) , where q'^ = (T"^(0) and = (0) (1 < « < A). Usually we write qQ instead 

of q^, and so we have the local chart (/3*'') ^ ([/) in T^Q with local coordinates (g^, g^, . . . , g^) . 
Local coordinates in T{T''Q) are denoted by [q^, qf^, . . . , g^; t^q^, ff^, . . . , u^) . 

Using these coordinates, the local expression of the canonical projections are p^. [q^, q^, . . . , q^) 
{qqi Qi-, ■ ■ ■ 1 Q^) ) a-nd then for the tangent maps T/oJ? : T{T''Q) — )• T{T^Q), we have the local ex- 
pression Tp^ {q^, , . . . , q^, v^, v^,..., v^) = [q^, qf,..., q^, v^, vf,..., v^) . 

If (t: M — )• Q is a curve in Q, the canonical lifting of a to T^Q is the curve ct*^: R — )• T^Q 
defined as a^it) = o"t'(0), where crt{s) = a{s + t), (that is, the A-jet lifting of a). If A = 1, we 
will write = a. 

Let V{p^_i) be the vertical sub-bundle of T^Q in T^~^Q. In the above coordinates, for 
every p G T^'Q and u G Vp{p'^_i), we have that its components are u = (O, . . . , 0, u^, . . . , v^^ . 
Furthermore, if ik-r+i- V{pr-i) T(T'^Q) is the canonical embedding, then 

ik^r+i {qo, ...,qk,Vr,---,Vk) = {<lo^ • • • , 9fc , 0, . . . , 0, t;;!^ , . . . , ) . 

Consider now the induced bundle of tjt-iq : T{T^~^Q) — J- T^^^Q by the canonical projection 
Pr-i' denoted by T'^Q x-pr-ig T(T''^^Q), which is a vector bundle over T^Q. We have the 
following commutative diagrams 

T'^Q XTr-iQT(T'-iQ)- - ^T(T'^-iQ) , T{T^Q) ^ T(T^-iQ) 
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Then, there exists a unique bundle morphism s^.^+i : T{T^Q) 
that the following diagram is commutative: 



T(T'=Q) 



T^Q x^r-iQ T(T^-iQ) such 




It is defined by Sfe_r+i('u) = [rrj^kqiu) ,T p'^_i{u)) , for every u G T(T'^Q). Its local expression is 

/4 A A A\ / A A A A A A \ 

Sk-r+1 [qo,---,qk,Vo,---,Vk) = [qo,..., qr-i,qr ,---,qk^VQ,..., V^_i} . 

As is a surjective map and Im (ik-r+i) = ker we have the exact sequence 

V{p^r-i) T(T'=Q) T'^Q XT.-1Q T(T'^-iQ) , 

which is called the (fc — r + 1) -fundamental exact sequence. In local coordinates, it is given by 

0^{q^,...,qi,vf,...,vi) "i^' {q^,...,qiA...Avf,...,vt) 

{q^,...,q^,v^,...,v^) 'tit' {q^, . . . ,q^;q^, . . . ,q^_i,v^, . . . ,v^_^) ^0 . 

Thus, we have A; exact sequences 



1st: ^ V{pt_i) ^ T{T'Q) ^ T'^Q Xt.-iq T(T'=-^Q) 



rth: ^ V^(Pfe-r) ^ T(T'=Q) ^ T^Q Xt^-.q T(T'=-'^Q) 



kth: ^ l^C^'^) ^ T(T'^Q) ^ T'=Q xq TQ ^ , 

where V{/3'') = ^(po) denotes the vertical subbundle of T'^Q on Q. These sequences can be 
connected by means of the connecting maps 



locally defined as 

,A „A „.A 



hk-r+i [qt ■■■^qt^t---^ ^tr) = (^g^, . . . , , 0, . . . , 0, ^v^, ^^^^^p^v^ 



k\ 



(k-r) 



These maps arc globally well-defined and arc vector bundle isomorphisms over T'^Q. Then we 
have the following connection between exact sequences: 







• Vip1_^) ^ T(T'=Q) ^ Ti^Q XT.-.Q T(T'=-'^Q) . 
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2.1.2 Higher-order canonical vector fields. Vertical endomorphisms and almost- 
tangent structures 



The canonical injection is the map 

cr^(O) 



where 



In local coordinates 



7: 



T(T'-iQ) 
7(0) 



, (1 < r < fc) , 



jr {Qo, • • • , C) = (C> • • • > it, Qt---, Qr) 



(1) 



(2) 



Then, the following composition allows us to define a vector field G 3i{T^Q), 

Idxjfe_^_n _,, _ ,, h 



T(T'=Q) ; 




that is, Ar = ik-r+i ° hk-r+i ° (Id x jk-r+i)- Prom the local expressions of ifc-r+i, ^fc-r+i and 

{k-ry.1k-r+ll' 



jk-r+1 we obtain that (g^, ...,q^) = (q^, . . . , g^, 0, . . . , 0, r! , (r + 1)! g^, . . . ^^-^ 



or what is equivalent 

k-r 



A. = ^ 

i=0 

In particular 



(k-r)] 



i=l i=0 



dgf+l 9q^ 9q^ 



d 



d 



Qk-r+l 



d 



_d_ 



Definition 1 The vector field A^ is the rth-canonical vector field in T^Q. In particular, Ai is 
called the Liouville vector field in T^Q. 



Remember that, if is a (fc + l)n-dimensional manifold, an almost-tangent structure of 
order A; in AT is an endomorphism J in TAT such that j'^'^^ = and rank J = kn. Then, T'^Q is 
endowed with a canonical almost-tangent structure. In fact: 



Definition 2 For 1 < r < k, let ik-r+i, hk-r+i, Sr be the morphisms of the fundamental exact 
sequences introduced above. The map 

Jr = ik-r+1 o hk-r+1 o : T(T'=g) T(T^Q) 

defined by the composition 



T(T'=g) 




T(T'=Q) 



is called the rth- vertical endomorphism ofT{T^Q). 
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Prom the local expressions of Sr, hk-r+i, ik-r+i we obtain that 





...,v^) = (q^,...,qtO,...Ar\v^,{r + l)\vt■.■ 



{k-r)\ 



kl 





As a consequence, the Ist-vertical endomorphism Ji defines an almost-tangent structure of 
order k in T^Q, which is called the canonical almost-tangent structure of T^Q. Then, any other 
vertical endomorphism Jr is obtained by composing Ji with itself r times. Furthermore, we 
have the following relation: 



As a consequence, starting from the Liouville vector field and the vertical endomorphisms, we 

can recover all the canonical vector fields. However, as all the vertical endomorphisms are 
obtained from Ji, we conclude that all the canonical structures in T'^Q are obtained from the 
Liouville vector field and the canonical almost-tangent structure. 

Consider now the dual maps J* of J^, 1 ^ r ^ k; that is, the endomorphisms in T*(T^(5), 
and their natural extensions to the exterior algebra /\(T*(T*^(5)) (also denoted by J*). Their 
action on the set of differential forms is given by 



for CO e i7P(T*^Q) and Xi, . . . ,Xp € X(T^Q), and for every / G C~(T^Q) wc write J*{f) = f. 
The endomorphism J* : f2{T'^Q) -> f2{T'^Q), 1 < r < k, is called the rth-vertical operator, and 
it is locally given by 



2.1.3 Vertical derivations and differentials. Tulczyjew's derivation 

The inner contraction of the vertical endomorphisms Jr with any differential p-form ui G 0P(T^Q) 
is the p-form i{Jr)oj defined as follows: for every Xi, . . . ,Xp e X{T^Q) 




iir + s^k-\-l 
ifr-t-s<fe-|-l 



j;uj{Xi, ...,Xp)= Uj{Jr{Xi), Jr{Xp)) , 



j:{f) = f 



for every / G C~(T'=Q) 




p 



i{JrMXi,...,Xp) = J2^{X: 




1=1 



and taking i{Jr)f = 0, for every / G C°°{T''Q), we can state: 
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Definition 3 The map 

UJ I > i{Jr)uJ 

is a derivation of degree in QiT^Q), which is called the rth-vertical derivation in QiT^Q). 
Its local expression is 

i{jrmt) 



0, if 2 < r 



Definition 4 The operator dj^ = [i(Jr),d] is a skew- derivation of degree 1, which is called the 
rth-vertical differential. 



Its local expression is given by 

k 



djrif) = Y.7—T^TrA^'ltr , for every /GC~(Tfcg) 
For 1 < r, s < k, we have that d j^d = —ddj^. 

In the set efe>o^?(T''Q), we can define the following equivalence relation: for co G Q{T''Q) 
and A E 0{T'''Q), 



Then we consider the quotient set f2 = fiiT^Q)/ ~, which is a commutative graded algebra. 

k^o 

In this set we can define the so-called Tulczyjew's derivation IlQl[l7], denoted by dx, as follows: 
for every / e C°°(T'=Q) we construct the function drf G C°^{T^+^Q) given by 

{dTf){a^+'m = (d^.(o)/)(jfc+i(<^''+^(0))) 

where jk+i '■ T'^+^Q ~^ T(T'^Q) is the canonical injection introduced in the Section [2.1.21 From 
the coordinate expression for jk+i, we obtain that 

dxf {qt • • • , it+i) = itijriilo^ ■■■,Qk)- 



This map dx extends to a derivation of degree in i? and, as d^d = ddx, it is determined by 
its action on functions and by the property dxidgf) = dq^^^. 

Furthermore, the maps i{Js), dj^, i{^s) and L(As) extend to 1? in a natural way. 



2.1.4 Higher-order semisprays 

Definition 5 A vector field X G X{T^Q) is a semispray of type r, I < r < k, if for every 



integral curve a of X, we have that, if j = (3'^ o a , then = p'^_j.^i o a (where j'' is 
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the canonical lifting of j to T'^ '^'^^Q). 




In particular, X G X{T'^Q) is a semispray of type 1 if for every integral curve a of X, we have 
that J = P'^ o a, then 7'^ = a. 



The local expression of a semispray of type r is 



d 



A d 



+ ■ ■ ■ + Qk-r+l-^-A ^ -^k-r+l 



dqtr 



Proposition 1 The following assertions are equivalent: 

1. A vector field X G X{T^Q) is a semispray of type r. 

2. Tp^_y, o X = that is, the following diagram commutes 

T(T*^Q) 

X 

3. Jr{X) = Ar. 



+ 




Obviously, every semispray of type r is a semispray of type s, for s > r. 

If X G X{T^Q) is a semispray of type r, a curve cr in Q is said to be a path or solution of 
X if cr'^ is an integral curve of X; that is, a'^ = X o a'^ , where a'^ denotes the canonical lifting 
of from T^Q to T(T'^Q). Then, in coordinates, a verifies the following system of differential 
equations of order A; + 1: 



d'^-^'+^a^ 



^k-^+^y^ dt''"' dt^ J 



da 



d''a\ 



^k+i^A 
dtk+^ 



Xt\o 



da 



d^c 



dt''"' dt^ 



Observe that, taking k = 1, then r = 1 and = Wtq, we recover the definition of the 

holonomic vector field (sode in TQ). So, semisprays of type 1 in T^Q are the analogue to the 
holonomic vector fields in TQ; that is, they are the vector fields whose integral curves are the 
canonical liftings to T^Q of curves on the basis Q. Their local expressions are 



X 



A 9 A d 
dqo dqf 



A 9 



9Qk-i 



^rA 9 



9<lt 



P.D. Prieto-Martmez, N. Roman-Roy: UniGed formalism for higher-order systems. 



10 



2.2 Lagrangian formalism 

Let Q be a n-dimensional differentiable manifold and £. G C°°{T''Q). We say that £ is a 
Lagrangian function of order k. 

Definition 6 The Lagrangian 1-form 6c G Q^{T'^^~^Q), associated to C is defined as 

k 

r! 



r=l 
')2frp2k-l/ 



Then, the Lagrangian 2- form, ojc G f2^(T^'' Q), associated to L is 

k 



= -dOc = Y.i-'^y^dJf'ddj^C 



r=l 



Observe that the Lagrangian 1-form is a semibasic form of type k in T^*^ . 

We assume that ojc has constant rank (we refer to this fact by saying that £ is a geometrically- 
admissible Lagrangian). 

Definition 7 The Lagrangian energy, Ec G C'^{T'^^~^Q), associated to C is defined as 

It is usual to write C instead of {p^''~^)*jC, and we will do this in the sequel. 
The coordinate expressions of these elements are 



r=l i=0 \^yr+i 
k k—r 



UC 



5:B-ir^^^Td(^lAd,-, 



r=l i=0 \^'ir+i 

k k—r 



Definition 8 A Lagrangian function C € C°°{T^Q) is said to be regular if uc is a symplectic 
form. Otherwise C is a singular Lagrangian. 

To say that £ is a regular Lagrangian is locally equivalent to saying that the Hessian matrix 
is regular at every point of T'^Q. 



Definition 9 A Lagrangian system of order k is a couple {T'^''~^Q, C), where Q represents the 
configuration space and L G C^i^^Q) is the Lagrangian function. It is said to he a regular 
(resp. singular j Lagrangian system if the Lagrangian function C is regular (resp. singular). 
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Thus, in the Lagrangian formalism, T^^~^Q represents the phase space of the system. The 
dynamical trajectories of the system are the integral curves of any vector field Xc G jC(T^'^~^(5) 
satisfying that: 

1. It is a solution to the equation 

i{Xc)u;c = dEc (5) 

2. It is a semispray of type 1 in T'^^^^Q. 

Equation ^ is the higher-order Lagrangian equation, and a vector field Xc solution to ([5]) (if 
it exists) is called a Lagrangian vector held of order k. If, in addition, X/^ satisfies condition 2, 
then it is called an Euler-Lagrange vector Geld of order k, and its integral curves on the base 
are solutions to the higher-order Euler-Lagrange equations. 

In natural coordinates of T^^~^Q, if 

2/c ]_ O ^1 



as 



k k T ^ k V k / c\2i /* \ ^ 



from ([5]) we obtain 



dq^dq^ 



2 r 2fc-3 



(-1)'= (/f.-r - ^ + Ec-ir'iT f S - E 



B B \ 



dq^dq^ ^ \dqi , .^^ 

(6) 

where the terms in brackets ( ) contain relations involving partial derivatives of the La- 
grangian and applications of dx, which for simplicity are not written. These are the local 
expressions of the Lagrangian equations for Xc. 

Now, if o": M — 7- T'^'^^^Q is an integral curve of Xc, from ([5]) we obtain that a must satisfy 
the Euler-Lagrange equation 

i{a){ijJc o a) = dEc o cr , (7) 

where a denotes the canonical lifting of a to T(T^^~^Q); and as Xc is a semispray of type 1, 
we have that a is the canonical lifting of a curve 7: M — )■ Q to T'^'^~^Q; that is, a = 

Now, if £ G C°°{T^Q) is a regular Lagrangian, then ujc is a symplectic form in T'^^'^Q, and 
as a consequence we have that: 

Theorem 1 Let {T'^^^^Q, C) be a regular Lagrangian system of order k. 



P.D. Prieto-Martinez, N. Roman-Roy: Unified formalism for higher-order systems. 



12 



1. There exists a unique Xc G X(T^'^ ^Q) which is a solution to the Lagrangian equation ^ 
and is a semispray of type 1 in T'^''~^Q. 

2. If J : ^ Q is an integral curve of Xc then a = j'^^^^ is a solution to the Euler-Lagrange 
equations; 

-dtw°' ^dr^w°' + 

If £ G C°°(T^(5) is a singular Lagrangian, then loc is a presymplectic form, so the existence 
and uniqueness of solutions to the Lagrangian equation ([5]) is not assured, except in special 
cases (for instance, when uc is a presymplectic horizontal structure [17J). In general, in the 
most favourable cases, equation ([5]) has solutions Xc G X(T^^~^Q) in some submanifold 5/ ^ 
T'^^~^Q, for which these vector fields solution are tangent. This submanifold is obtained by 
applying the well-known constraint algorithms (see, for instance, [22\ [2T\ I30j). Nevertheless, 
these vector fields solution are not necessarily semisprays of type 1 on Sf, but only on the 
points of another submanifold Mj Sf ^ T'^^~^Q (see |21i I30j). On the points of this 
last submanifold, the integral curves of Xc G X{T^^~^Q) are solutions to the Euler-Lagrange 
equations dS]). 

A detailed study of higher-order singular Lagrangian systems can be found in [231 124j . 



2.3 Hamiltonian formalism 

Definition 10 Let (T'^^^^Q, C) be a Lagrangian system. The Legendre-Ostrogradsky map (or 
generalized Legendre mapj associated to C is the map TL: T^'^^^Q — t- T*(T'^^^Q) defined as 
follows: for every u G T(T'^^~^Q), 



ec{u) = (Tpf_-i(n),-F£(rT2.-iQ(u))) 



This map verifies that -Kr^k-iqoFC = p^_i , where vTrpfc^ig : T*(T^' '^Q) — )■ T'^' is the nat- 
ural projection. Furthermore, if 6'fc_i G Q^{T*{T^~^Q) and oj^-i = -dOk-i G Q"^ {T* {T''-^ Q)) 
are the canonical 1 and 2 forms of the cotangent bundle T*(T^-iQ), we have that 

FCOk-i = 0c , FC*uJk~i = oJc ■ 
Given a local natural chart in T'^^~^Q, we can define the following local functions 

k—r 

Pa' = E(-l)'^^ 

Observe that 



^ dL ^ 



and hence 

Pa~' = §^ - Mpa) , l<r<A;-l. (9) 
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Thus, bearing in mind the local expression Q of the form 9c, we can write 6c = X]r=i P^A ^'^Qt^-i^ 
and we obtain that the expression in natural coordinates of the map J^C is 

J^C {qt qt,---, Q2k-l) = {qo, qt,---, qtl,PA,PA, ■ ■ ■ ,Pa~') ' with p\oTC= p\ . 

C is a. regular Lagrangian if, and only if, FC: T^^^^Q — )• T*{T^^^Q) is a local diffeomor- 
phism. 

As a consequence of this, we have that, if £ is a regular Lagrangian, then the set 
< i < k — 1, is a set of local coordinates in T^'^^^Q, and (p^) are called the Jacobi-Ostrogradsky 
momentum coordinates. 

Observe that the relation ([9]) means that we can recover all the Jacobi-Ostrogadsky momen- 
tum coordinates from the set (j5^~^). 

Definition 11 £ G C°°{T''Q) is said to be a hyperregular Lagrangian of order k ij TL is a 
global dijjeomorphism. Then, {T'^^~^Q, C) is a hyperregular Lagrangian system of order k. 

As TTrpfc-iQ o = plf^Ti^, this condition is equivalent to demanding that the restriction of 
^2k~i . rp2fc-ig _^ t'^-^Q to every fibre be one-to-one. 

In order to explain the construction of the canonical Hamiltonian formalism of a Lagrangian 
higher-order system, we first consider the case of hyperregular systems (the regular case is the 
same, but restricting on the suitable open submanifolds where J-C is a local diffeomorphism) . 

So, {T'^^~^Q, C) being a hyperregular Lagrangian system, since FC is a diffeomorphism, there 
exists a unique function h e C°°(T*(T'=-1Q)) such that FC*h = Ec, which is called the Hamil- 
tonian function associated to this system. Then the triad {T*{T'^~^Q),uJk-i, h) is called the 
canonical Hamiltonian system associated to the hyperregular Lagrangian system [T^^~^Q,C). 
Thus, in the Hamiltonian formalism, T*(T^~^Q) represents the phase space of the system. 

The dynamical trajectories of the system are the integral curves of a vector field Xh € 
X{T*{T^~^Q)) which is a solution to the Hamilton equation 

iiXh)u;k-i = dh . (10) 

As ujk-i is symplectic, there is a unique vector field solution to this equation, and it is called 
the Hamiltonian vector field. 

In natural coordinates of T*{T^^^Q), {qf,p\) (with 0<i<k — 1; l<A<n), taking 
d ^ ^ d dh ^ A dh 

that 

fA_ ^ i ^ dh 

dp\ ' dqf • 

Now, if (t: M — )• T*{T'^^^Q) is an integral curve of Xh, we have that a must satisfy the 
Hamiltonian equation 

i((7)(a;fc_i o a) = dhoa , 

and, if a{t) = {qf{t),p^^{t)) in coordinates, it gives the classical expression of the Hamilton 
equations: 

dqf- dh dp\ dh 

= — ^oa , = T°cr- 

dt dp\ dt dqf 



= fi^Tr^+9ATr-r^ = TT^f^qf + tt-t^Pa^ and ujk-i = dqfAdp\, from ([10]) we obtain 

dqf dp\ dqf dp\ 



P.D. Prieto-Martinez, N. Roman-Roy: Unified formalism for higher-order systems. 



14 



For the case of singular higher-order Lagrangian systems, in general there is no way to asso- 
ciate a canonical Hamiltonian formalism, unless some minimal regularity condition are imposed 
[23]. In particular: 

Definition 12 A Lagrangian C € C°°{T^Q) is said to be an almost-regular Lagrangian function 
of order k if: 

1. TOT'^^-^Q) = Po is a closed suhmanifold of T* {T''-^ Q) . 
(We denote the natural embedding by jp^ : Pq ^ T*(T^^^Q)). 

2. T C is a surjective submersion on its image. 

3. For every p G T'^^~^Q, the fibers J- C^^ {J- C{p)) are connected submanifolds of T^^^^Q . 
Then {T'^^^^Q,C) is an almost-regular Lagrangian system of order k. 

Denoting the map defined by = jp^ o TCo by TCo- T'^^^^Q — Pq, we have that the 
Lagrangian energy Ec is a .FjCo-projectable function, and then there is a unique function ho € 
C~(Po) such that J^Qho = El (see [23j). 

This ho is the canonical Hamiltonian function of the almost-regular Lagrangian system and, 
taking lOq = jp^Wfc-ii the triad {Po,uJo, ho) is the canonical Hamiltonian system associated to the 
almost regular Lagrangian system (T'^^~^Q, C). For this system we have the Hamilton equation 

iiXh^)ojo = dho , Xh^ e X(Po) . (11) 

As ojo is, in general, a presymplectic form, in the best cases, this equation has some vector 
field solution only on the points of some submanifold Pj ^ Po ^ T*(T'^~^Q), for which 
Xh^ is tangent to Pf. This vector field is not unique, in general. It can be proved that Pf = 
J-C{Sf), where Sf T^^~^Q is the submanifold where there are vector field solutions to 
the Lagrangian equation ([5]) which are tangent to Sf (see the above section). Furthermore, 
as J- Co is a submersion, for every vector field X^^ € X(T*(T'^~^Q)) which is a solution to 
the Hamilton equation (fTT]l on Pf, and tangent to Pf, there exists some semispray of type 1, 
Xc G X{T'^'^^^Q), which is a solution of the Euler-Lagrange equation on Sf, and tangent to Sf, 
such that TCo*Xc = X^^. This J^/^o-projectable semispray of type 1 could be defined only on 
the points of another submanifold Mf ■— >• Sf. (See [231 [2l] for a detailed exposition of all these 
topics) . 

3 Skinner-Rusk unified formalism 

3.1 Unified phase space. Geometric and dynamical structures 

Let C € C°°(T^Q) be the Lagrangian function of order k of the system. First we construct the 
unified phase space 

(the fiber product of the above bundles), which is endowed with the canonical projections 

pri : T2'^-iQxTfc-iQT*(T^"iQ) ^ T^'^-Iq ; prg : T^''~^Qx^k-iQT* (T'^'^Q) ^ T*(T'="1Q) , 
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and also with the canonical projections onto T'^ ^Q. So we have the diagram: 





If {U, Qq) is a local chart of coordinates in Q, denoting by {{13'^'^^^)^^ (U); q^, Qi, . . . , 
and ((7rrpfc_iQO^'^~i)~i([7);g^,gfj^, . . . ,q^_i,p%p% . . . the induced charts in T^'^^^Q and 

in T*{T'^~^Q), respectively, we have that (qq, . . . , q^_i] q^, ■ ■ ■ , 'l2k-i^P\^ ■ ■ ■ ^P^A^) ^'^^ 
ural coordinates in the suitable open domain W C W. Note that dim(W) = 2>kn. 

The bundle W is endowed with some canonical geometric structures. First, let E 
J72(T*(T^-^Q)) be the canonical symplectic form oiT*{T:^-^Q). Then we define 

which is a presymplectic form in W, whose local expression is 

= pr^ Wfc-i = pr^ {dqt A dp^) = dg,^ A dp\ . (12) 

Observe that 

kerJ^ = (|^,...,^)=^^(-^Hw). (13) 
The second relevant canonical structure in W is the following: 



Definition 13 Let p S T^'^ ^Q, its projection q = p^_i{p) to ^Q, and a covector Uq E 
T*(T'=-iQ). The coupling function C G C°°(>V) is defined as follows: 

C:T2^-1Qxt.-iqT*(T'^-1Q) R 

(P,ag) ' — > {aq\jk{p)q) ' 

where j^: T^^~^Q — )■ T(T'^~^(5) is i/ze canonical injection introduced in ([I]), jk{p)q is the cor- 
responding tangent vector to T^^^Q in q, and {uq \ jk{p)q) = Ciq{jk{p)q) denotes the canonical 
pairing between vectors of Tq(T^^^Q) and covectors of T*{T^^^Q). 



Note that, in this case, jk '■ T Q — t- T(T Q) is a diffeomorphism 



In local coordinates, if p = {q^, q^_^, g^, • • • , g^-i)' t^^^n q = pf_^{p) = {q^, qk_i} 
and bearing in mind the local expression ([2]) of jk, we have jfc(p) = (^o^, • • • , Qt-ii^i^ • • • ' ^k) 

d 

Therefore if jk{p)q 



Qf+1 



G Tq(T Q), and if aq = p\ dq^ we obtain the following 



local expression for the coupling function C 



C{p, aq) = {aq I jk{p)q) = ( dqt 



1^ 

li+l 



dqf 



(15) 
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Observe that, if A; = 1, the map ji : TQ — )■ TQ is the identity on TQ, and we recover the 
standard canonical couphng between vectors in TpQ and covectors in T*Q. 

Using the couphng function, given a Lagrangian function C G C°°{T''Q), we can define the 
Hamiltonian function H G C°°(W) as 

H = C-{pf~^o^,,YC, (16) 

whose coordinate expression is 

H=fj,qt,^-C{qt....qt) • (17) 
Now, (W, Vt, H) is a presymplectic Hamiltonian system. 

Finally, in order to give a complete description of the dynamics of higher-order Lagrangian 
systems, we need to introduce the following concept: 



Definition 14 A vector field X G X(W) is said to be a semispray of type rinW if, for every 
integral curve a: I C M — ?■ W of X, the curve ai = prj^ofj: / — t- T'^^~^Q satisfies that, if 



In particular, X G X(>V) is a semispray of type 1 if 72^^-1 =(71. 



The local expression of a semispray of type r in W is 

2 A' 1 c\ 2 1 1 

■1=0 ^* i=2k—r * i=0 ^ 

and, in particular, for a semispray of type 1 in W we have 

2/c 2 r\ k 1 c-~. 

^ 2^ ^i+l + "^2fc-lT— 4 1- / .^A^TT ■ 



3.2 Dynamical vector fields 

3.2.1 Dynamics in W = T^'^-iQ x^k-iQ T*(T^-iQ) 

As we know, the dynamical equation of the presymplectic Hamiltonian system (W, H) is 
geometrically written as 

i{X)n = dH ; X e X(W) . (18) 

Then, according to |22] we have: 



Proposition 2 Given the presymplectic Hamiltonian system {W,il.,H), a solution X G X{W) 
to equation ^8\} exists only on the points of the submanifold Wo ^ W defined by 

Wo = {peW: ^{p) = {i{Y)dH){p) = , Vy G ker f)} . (19) 



We denote by jo ■ Wo ^ W the natural embedding and by (W) the set of vector fields 
in W at support on Wo- We have the following characterization of Wo- 
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Proposition 3 The submanifold Wo ^ W is the graph of the Legendre-Ostrogradsky map de- 
fined by C; that is, Wo = graph 



(Proof) As Wo is defined by it suffices to prove that the constraints defining Wo are those 
defining the graph of the Legendre-Ostrogradsky map associated to C. We make this calculation 
in coordinates. Taking the local expression (jl7p of the Hamiltonian function H E C°°(W), we 
have 

h 1 k 

dH = ^{qf^,dp\+p\dqt,^) - ^ -^dqf , 

i=0 1=0 



and using the local basis of ker Q given in (jl3p . we obtain that the equations defining the 
submanifold Wo are 



i(Y)dH = 0^p'X-'--^ = 



Observe that these expressions relate the momentum coordinates p^ ^ with the Jacobi-Ostrogadsky 
dC 

functions p^a~^ = — j, and so we obtain the last group of equations of the Legendre-Ostrogradsky 

map. Furthermore, in Section 12.31 we have seen that the other Jacobi-Ostrogradsky functions 
P^A^ (l<r<fe — 1) satisfy the relations ([9]). Thus we can consider that Wo is the graph of a 
map 

F: T^'^-iQ — > T*{T^^^Q) 

(qf) ^ {q^,...,qt„P%---,p'A-') 
which we identify with the Legendre-Ostrogradsky map by making the identification = P^^^ ■ 



Hence, we look for vector fields X £ Xw^ (W) which are solutions to equation (fT8|) at support 
on Wo] that is 

[i{X)n-dH]\^^=0. (20) 
In natural coordinates a generic vector field in X{W) is 



i=0 



bearing in mind the local expressions of $7 and dH, from (jlSp . we obtain the following system 
of {2k + l)n equations 



= iti , (21) 



GA = TrA ' G\ = -^-p'X=dT{PA), (22) 
oqQ dqf 



dC 



Pa-^ = ^^ (23) 
dq^ 

where ^ i ^ A; - 1 in (gl]) and 1 ^ ? s; - 1 in Therefore 

k — ]_ r\ '2,k 1 r\ r\ r\ k — 1 r\ 



We can observe that equations (j23p are just a compatibility condition that, together with the 
other conditions for the momenta, say that the vector fields X exist only with support on 
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the submanifold defined by the graph of the Legendre-Ostrogradsky map. So we recover, in 
coordinates, the result stated in Propositions [2] and [3l Furthermore, this local expression shows 
that X is a semispray of type k'mW. 



The component functions Ff-, k ^ i ^ 2k — 1, are undetermined. Nevertheless, we must 
study the tangency of X to the submanifold Wo', that is, we have to impose that L(X)^|yy = 
^(C)lw ~ ^' ^'-'^ every constraint function ^ defining Wo- So, taking into account Prop. [3l 
these conditions lead to 



and, from here, we obtain the following kn equations 
B B 




dq 



(T-iB B \ ^ ^ ( i~iB B \ 1 { d ^ 



^1 £^ A;— 3 

- (ik-\) Br. A - 5Z (^-fc+i ~ 9fc+i+l) ( ) = 

where the terms in brackets ( ) contain relations involving partial derivatives of the La- 

grangian and applications of which for simplicity are not written. These are just the La- 
grangian equations for the components of X, as we have seen in ([6]). These equations can be 
compatible or not, and a sufficient condition to ensure compatibility is the regularity of the 
Lagrangian function. In particular, we have: 

Proposition 4 // £ G C°°(T'^(5) is a regular Lagrangian function, then there exists a unique 
vector field X G X>v„(W) which is a solution to equation i20\) : it is tangent to Wo, and is a 
semispray of type linW. 

(Proof) As the Lagrangian function C is regular, the Hessian matrix I „ . I is regular at 

ydqj^dq^J 

every point, and this allows us to solve the above k systems of n equations (I25p determining all 
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the functions F^^ uniquely, as follows 

F-*=qf+i , (fc<i<2t-2) (26) 

(-1)' (Fi., - ,r feVO) ^ . h-^'^ (^) ^ . 

In this way, the tangency condition holds for X at every point on Wo- Furthermore, the equalities 
(f26l) show that X is a semispray of type 1 in W ■ 



However, if C is not regular, the equations (j25p can be compatible or not. In the most 
favourable cases, there is a submanifold W/ ^ Wo (it could be W/ = Wo) such that there exist 
vector fields X G XyVo(W), tangent to W/, which are solutions to the equation 

[i{X)n-dH]\y,^=Q . (27) 

In this case, the equations ()25p are not compatible, and the compatibility condition gives rise to 
new constraints. 



3.2.2 Dynamics in T^*^ 

Now we study how to recover the Lagrangian dynamics from the dynamics in the unified for- 
malism, using the dynamical vector fields. 

First we have the following results: 

Proposition 5 The map pr;^ = pr;^ ojo : is a diffeomorphism. 

(Proof) As Wo = graph. J^C, we have that T'^''~^Q ~ Wo- Furthermore, pfi is a surjective 
submersion and, by the equality between dimensions, it is also an injective immersion and hence 
it is a diffeomorphism. ■ 

Lemma 1 If uJk-i £ f2'^{T*{T^~^Q)) is the canonical symplectic 2-form in T*(T^~^(5), and 
cj£ = J^C*UJk-i is the Lagrangian 2-form, then VL = prj ujc- 

[Proof) In fact, 

pr* ujc = pT\{TC*uJk-i) = {TC o pri)*u;fc_i = pr2 LOk-i = . 



Lemma 2 There exists a unique function Ec € Qoo(^^'ik Iq-^ such that pr^ E'^ = H. This 
function Ec is the Lagrangian energy. 

(Proof) As pr^ is a diffeomorphism, we can define the function E^ = (pr]~^ o jo)* LI G 
C°°(T^'^^^(5), which obviously verifies that prj Ec = LI. 

In order to prove that Ec is the Lagrangian energy defined previously, we calculate its local 
expression in coordinates. Thus, from (I17p we obtain that 

fc-i 



WlEc = H = Y^PAlt+i - ^Qo, ■■■,Qk 



i=0 
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but Wo ^ W is the graph of the Legendre-Ostrogradsky map, and by Prop. [3] we have p\ 



i=o \^lUi+j, 



k 1 k — X — 1 / r\ /-I 

OL 



i=0 i=o V^^i+i+i/ 

= E E <it{-^y<^T {^] - ^iit ■■■^it)- 

i=l i=0 K^'^i+jJ 
Now, as pfi = pr^ ojg and pr* qf = qf, we obtain finahy 

= E E ^fi-^yd^ ii^] - ^(90^ ■ ■ ■ , Q^) 

which is the local expression (j3D of the Lagrangian energy. ■ 

Using these results, we can recover an Euler-Lagrange vector field in T'^^~^Q starting from 
a vector field X G Xw;q(W) tangent to Wo, a solution to (|20l) . First we have: 

Lemma 3 Let X £ X(W) be a vector field tangent to Wo- Then there exists a unique vector 
field Xc G X{T'^^~^Q) such that Xc o pi^ oj^ = Tpr;^ oX o j^. 

(Proof) As X G X{W) is tangent to Wo, there exists a vector field Xq G 3l{Wo) such that 
Tjo o Xo = X o jo- Furthermore, as pr^ is a diffeomorphism, there is a unique vector field 
Xc G Xir'^^^^Q) which is pr;^-related with Xo] that is, Xc o pf;^ = Tpf]^ o Xo- Then 

Xc o pri ojo = Xcopr[ = Tpf 1 oXo = T pii oT jo o Xo = T pr^ oX o 



And as a consequence we obtain: 



Theorem 2 Let X G X>v„ ( W) be a vector field solution to equation (20\) and tangent to Wo ( at 
least on the points of a submanifold W/ ^ Wo)- Then there exists a unique semispray of type 
k, Xc G X{T'^^^^Q), which is a solution to the equation 

i{Xc)oJc -dEc = (28) 

(at least on the points of Sf = pTi{Wf))- Ln addition, if C G C°°(T'^Q) is a regular Lagrangian, 
then Xc is a semispray of type 1, and hence it is the Euler-Lagrange vector field. 

Conversely, if Xc G X(T'^'^^^Q) is a semispray of type k (resp., of type 1), which is a solution to 
equation i28]) (at least on the points of a submanifold Sf T'^^'^^Q), then there exists a unique 
vector field X G XyVo(W) which is a solution to equation ![2U\) (at least on W/ = pfj~^(5/) 
Wo ^ W), and it is a semispray of type kinW (resp., of type 1 ). 

(Proof) Applying Lemmas [TJ [21 and El we have: 

0= [i{X)n - dH]\^^ = [i{X)pilujc-dpTlEc]\y^^=pTl [i{Xc)u;c-dEc]\^^ , 
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but, as pii is a surjective submersion, this is equivalent to 

= [i{Xc)u;c - di^/:]|p,^(w„) = [i{Xc)uJc - di^rllx^.-ig = , 

smce pr^ (Wo) = T^^-iQ. The converse is immediate, reversing this reasoning. 

In order to prove that Xc is a semispray of type k, we proceed in coordinates. From the 
local expression ()24p for the vector field X (where the functions are the solutions of the 
equations (|25|) ). and using Lemma [31 we obtain that the local expression of Xjr G X{T^^^^Q) is 



fc-i g 2k-i Q 

and then 

so Xc is a semispray of type k in T^^^^Q. 

Finally, if £ G G^iT^Q) is a regular Lagrangian, equations ()25p become (j26p . and hence the 
local expression of X is 



Therefore 



and then Ji{Xc) = (i + l)g^i ^ = Ai, which shows that is a semispray of type 1 
i=o 

in T^'^-iQ. ■ 

Remarks: 



• It is important to point out that, if C is not a regular Lagrangian, then X is a semispray 
of type k in W, but not necessarily a semispray of type 1. This means that Xc is a 
Lagrangian vector field, but it is not necessarily an Euler-Lagrange vector field (it is not a 
semispray of type 1, but just a semispray of type k). Thus, for singular Lagrangians, this 
must be imposed as an additional condition in order that the integral curves of Xc verify 
the Euler-Lagrange equations. This is a different from the case of first-order dynamical 
systems (A; = 1), where this condition {Xc is a semispray of type 1; that is, a holonomic 
vector field) is obtained straightforwardly in the unified formalism. 

In general, only in the most interesting cases have we assured the existence of a submanifold 
W/ ^ Wo and vector fields X G dLyj^iyV) tangent to W/ which are solutions to the 
equation (pT]) . Then, considering the submanifold Sf = pr^(>V/) T'^^~^Q, in the best 
cases (see [QESlIS^), we have that those Euler-Lagrange vector fields Xc exist, perhaps on 
another submanifold Mf ^ Sf where they are tangent, and are solutions to the equation 



[ixc^c-dEc]\M =0 



(29) 
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• Observe also that Theorem [5] states that there is a one-to-one correspondence between 
vector fields X G Xvv;„(>V) which are solutions to equation ([20]) and Xc € X(T^'^~^Q) 
solutions to ([281) . but not uniqueness, unless C is regular. In fact: 

Corollary 1 // £ G C'^{T^Q) is a regular Lagrangian, then there is a unique X G Xyv;^(W) 
tangent to Wo which is a solution to equation [20\). and it is a semispray of type 1. 

(Proof) As C is regular, by Proposition[T]there is a unique semispray oftypel,X£ G X{T^^-^Q) 
which is a solution to equation (|28p on T'^^^^Q. Then, by Theorem [21 there is a unique X G 
jCw^(W), tangent to Wo, which is a solution to ([20l) on Wo- ■ 

3.2.3 Dynamics in T*(T'=-^Q) 

Hyperregular and regular Lagrangians 

In order to recover the Hamiltonian formalism, we distinguish between the regular and non- 
regular cases. We start with the regular case, although by simplicity we analyze the hyperregular 
case (the regular case is recovered from this by restriction on the corresponding open sets where 
the Legendre-Ostrogadsky map is a local diffeomorphism) . For this case we have the following 
commutative diagram 



TW 




Theorem 3 Let C € C°°(T^Q) be a hyperregular Lagrangian, h G C°°(T*(T^-iQ)) the Hamil- 
tonian function such that TC*h = Ec, and X G X>v^(W) the vector field solution to the equation 
(2^) . tangent to Wo- Then, there exists a unique vector field X^ = TC^^Xn G X(T*(T'^^-'^(5)) 
which is a solution to the equation 



i{Xh)uJk-i -dh = 



(30) 
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Conversely, if G X{T*{T^ ^Q)) is a solution to equation /i^U\) . then there exists a unique 
vector field X G X>Vo(W), tangent to Wo, which is a solution to equation 120\) . 

(Proof) If C is hiperregular, then = J-C o pfj^ is a diffeomorphism, since it is a composition 
of diffeomorphisms; then there exists a unique vector field Xq G X(Wo) such that pr2*-^o = X^, 
and there is a unique X £ Xy\;^{W) such that jo^^Xo = X\y^^. 

Now, as FC*h = Ec, by applying Lemma [2] we have that ]ii\{FC*{h)) = pr| £'£ = H; but 
J-C o pr^ = pr2, and then pr2 h = H. Therefore, by the definition of fi, we have 

0= [iiX)n-dH]\^^= [^iX)prlu;k-l-dpr*^h]\^^=prl [^iXh)u;k-l - dh]\^^ . 
However, as pr2 is a surjective submersion and pr2 (Wo) = T*{T^^^Q), we finally obtain that 
= [i{Xh)uJk-i - cl/i]|pi.2(Wo) = [i{Xh)^k-i - d/i]|T*{Tfe-iQ) 



Singular (almost-regular) Lagrangians 

Remember that, for almost-regular Lagrangians, only in the most interesting cases have we 
assured the existence of a submanifold W/ ^ Wo and vector fields X G Xh'()(W) tangent to W/ 
which are solutions to equation ()27p . In this case, the dynamical vector fields in the Hamiltonian 
formalism cannot be obtained straightforwardly from the solutions in the unified formalism, but 
rather by passing through the Lagrangian formalism and using the Legendre-Ostrogadsky map. 

Thus, we can consider the submanifolds Sf = pri(W/) T^'^^^Q and Pf = pr2(W/) = 
TC{Sf) ^ T*(T'=-^Q). Then, using Theorem El from the vector fields X G 3ew„(W) we obtain 
the corresponding Xc G X(T^^~^Q), and from these the semisprays of type 1 (if they exist) 
which are perhaps defined on a submanifold Mf ^ Sf, are tangent to Mf and are solutions to 
equation (|29p . So we have the following commutative diagram 



W 




Now, it is proved ([24J) that there are Euler-Lagrange vector fields (perhaps only on the 
points of another submanifold Mf ^ Mf), which are J^i2-projectable on Pf = FC{Sf) ^ Pq ^ 
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T*{T''~^Q). These vector fields X^, = TL^Xc G X(T*(T'=-1Q)) are tangent to P/ and are 
solutions to the equation 

{^Xn,^o-^h^\p^ =0 . 

Conversely, as TLq is a submersion, for every solution Xh, G X(T*(T*'-iQ)) to the last 
equation, there is a semispray of type 1, X^ G jt(T^'^~^(5), such that FCq^Xc = Xh^, and we 
can recover solutions to equation (|27|) using Theorem [2J 

3.3 Integral curves 

After studying the vector fields which are solutions to the dynamical equations, we analyze their 
integral curves, showing how to recover the Lagrangian and Hamiltonian dynamical trajectories 
from the dynamical trajectories in the unified formalism. 

Let X G Xwo(W) be a vector field tangent to Wo which is a solution to equation (pOj) . and 
let fj : / C M — 7- W be an integral curve of X, on Wo- As it = Xoa, this means that the following 
equation holds 

i{a){VL o a) = dH o a . (31) 

Furthermore, if ctq : / — )• Wo is a curve on Wo such that joOf^o = a, we have that CTo is an integral 
curve of the vector field Xo G X(Wo) associated to X, and (Tq = Xq o Uo- 

In local coordinates, if a{t) = {qf{t),p'^^{t)), we have that 

qf{t)=qf+ioa {O^i^k-l) ; qf{t)=F,^oa {k ^ i ^ 2k - 1) 

dC 

where F/^ are solutions to equations ([25]) . 

Now, for the Lagrangian dynamical trajectories we have the following result: 

Proposition 6 Let cj : I C K — )• W be an integral curve of a vector field X solution to i20\). on 
Wo- Then the curve ac = pi^ oa: I ^ T^^~^Q integral curve of Xc- 

{Proof) As (T = jo o (To, using that Tjo o Xq = X o jo and that T pr;^ oX = Xc o pr;^, we have 

(T£ = pr7°S- = pri ojo o fjo = T pr^ oTjo o o-q = T pr^ oTjo oXooao 
= T pr^ oX o jo oao = Xco pr^ ojo o ao = X^o ac ■ 



Corollary 2 If L G C°°(T^Q) is a regular Lagrangian, then the curve ac = pri oo" : / — t- 

is the canonical lifting of a curve on Q; that is, there exists 7: / C M — )■ Q such that ac = "y'^^~^ ■ 

{Proof) It is a straighforward consequence of Proposition [6] and Theorem [2J ■ 
And for the Hamiltonian trajectories, we have: 

Proposition 7 Let a : / C M — s- W he an integral curve of a vector field X solution to i20\) . on 
Wo- Then the curve ah = TLo ac'- I ^ T*{T^~^Q) is an integral curve of X^ = FC^,{Xc)- 
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(Proof) Given that ac is an integral curve of Xc, Proposition El and the definitions of and 
(T/i, we obtain 

dh = J^C o ac = TTC o ac = TTC o Xc ° oc = X^ o J^C o cr£ = Xh o ah ■ 
Thus ah is an integral curve of Xh- ■ 
The relation among all these integral curves is summarized in the following diagram 




Remark: Observe that in Propositions [6] and [7] we make no assumption on the regularity of 
the system. The only considerations in the almost-regular case are that, in general, the curves 
are defined in some submanifolds which are determined by the constraint algorithm, and that ac 
is not necessarily the lifting of any curve in Q and this condition must be imposed. In particular: 

• If the Lagrangian is regular (or hiperregular), then Im(cj) C Wo, Im (cr^) C T^^^^Q and 
ImK) CT*(T'=-iQ). 

• If the Lagrangian is almost-regular, then Im (u) C W/ ^ Wo, lm.{ac) <Z Sf ^ T^^^^Q 
and Im(cj/,) C Pf ^ Po ^ T*{T''~^Q). 

4 Examples 

4.1 The Pais-Uhlenbeck oscillator 

The Pais-Uhlenbeck oscillator is one of the simplest (regular) systems that can be used to 
explore the features of higher order dynamical systems, and has been analyzed in detail in many 
publications [32^ |28] . Here we study it using the unified formalism. 

The configuration space for this system is a 1-dimensional smooth manifold Q with local 
coordinate {qo). Taking natural coordinates in the higher-order tangent bundles over Q, the 
second-order Lagrangian function C G C°°(T^(5) for this system is locally given by 



^qo,qi,q2) = ^ {ql- ^^'^ql- iql) 
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where 7 is some nonzero real constant, and a; is a real constant. £ is a regular Lagrangian 
function, since the Hessian matrix of C with respect to (72 is 



dq2dq2 



-7 



which has maximum rank, since we assume that 7 is nonzero. Notice that, if we take 7 = 0, 
then C becomes a first-order regular Lagrangian function, and thus it is a nonsense to study this 
system using the higher-order unified formalism. 

As this is a second-order dynamical system, the phase space that we consider is 

>V = T3Qxtq T*(TQ) 



T*(TQ) 




Denoting the canonical symplectic form by G f2'^{T*{TQ)), we define the presymplectic 
form Q G pr2iOi G /7^(W) with the local expression 

O = dgo A dp° + dqi A dp^ , 

The Hamiltonian function H G C°°(W) in the unified formalism is H = C — (p^ o pr^)*£, where 
C is the coupling function, whose local expression is 

C(go,gi,g2,g3,p°,p^) =P°Ql +P^92 • 
and then the Hamiltonian function can be written locally 

Hiqo,qi,q2,q3,P°,P^) =P°qi+P^q2 - ^ {qf -^"^ql-iql) 



As stated in the above sections, we can describe the dynamics for this system in terms of 
the integral curves of vector fields X G 3C(W) which are solutions to equation (jlSp . If we take a 
generic vector field X in W, given locally by 

d d d d d d 

X = /ot; 1- /l^ 1- ^2-^ 1- -^3^ 1- G -I- G — — j-, 

oqo oqi uq2 oq^ op" op^ 

taking into acount that 

dH = u?qodqQ + {p° - qi)dqi + {p^ + 792)dg2 + gidp° + q2dp^ , 

from the dynamical equation i{X)i} = dH, we obtain the following system of linear equations 
for the coefficients of X 

fo = qi (32) 

/i = 92 (33) 

= -uj\o (34) 

G' = qi-p'> (35) 

p^ + 7*72 = (36) 
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Equations ([55]) and give us the condition of semispray of type 2 for the vector field X. 
Furthermore, equation (j36p is an algebraic relation stating that the vector field X is defined 
along a submanifold Wo that can be identified with the graph of the Legendre-Ostrogradsky 
map, as we have seen in Propositions [2] and [3l Thus, using ([32]) . ([33]) . and ([35]) . the vector 
field X is given locally by 

X = qi—+ q2— +F2— + F3—- W^QOTH) + ill - P° • 37 
oqo oqi oq2 oqs op^ dpi 



As our goal is to recover the Lagrangian and Hamiltonian solutions from the vector field 
X, we must require X to be a semispray of type 1. Nevertheless, as £ is a regular Lagrangian 
function, this condition is naturally deduced from the formalism, as we have seen in (|25|) . 

Notice that the functions F2 and F^ in ()37p are not determined until the tangency of the 
vector field X on Wo is required. Recall that the Legendre-Ostrogradsky transformation is the 
map FC: T^Q — > T*(TQ) given in local coordinates by 

, (dC\ dC ^ , 

-7^2 

and, as 7 7^ 0, we see that £ is a regular Lagrangian since FC is a (local) diffeomorphism. Then, 
the submanifold Wo = graph FC is defined by 

Wo = {pG W:eo(p)=6(p) = 0} , 

where = — FC*p^, r = 1,2. The diagram for this situation is 



FC*p' = ^ 
oqi 

dC 



FC*p' = — 
oq2 




Next we compute the tangency condition for X £ X(W) given locally by (f37l) on the sub- 
manifold Wo ^ W, by checking if the following identities hold 

L(X)eo|w„=0 , L(X)ei|wo=0. 

As we have seen in Section [3.2. 11 these equations give us the Lagrangian equations for the vector 
field X; that is, on the points of Wo we obtain 

L(X)eo = -oj\o - 92 - 7^3 = (38) 
L(X)ei = 7 (F2 - gs) = . (39) 

Equation (j39p gives us the condition of semispray of type 1 for the vector field X (recall that 
7/0), and equation ([38]) is the Euler-Lagrange equation for the vector field X. Notice that, as 
7 is nonzero, these equations give us a unique solution for F2 and F3. Thus, there is a unique 
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vector field X £ X(W) solution to the equation — di7]|j^^ = which is tangent to the 

submanifold Wo ^ W, and it is given locally by 

d d d I , o \ d 2 9 , ()\ 9 

oqo oqi oq2 7 093 op'^ opi 

Then, if cr : M —>■ W is an integral curve of X locally given by 

a{t) = {qo{t),qi{t),q2{t),q3{t),p'^{t),p\t)) , (40) 
and its component functions are solutions to the system 

(41) 
(42) 
(43) 

uj\o{t) + q2{t)) ; (44) 

lo{t); (45) 
-/(t). (46) 

Finally we recover the Lagrangian and Hamiltonian solutions for this system. For the La- 
grangian solutions, as we have shown in Lemma [3] and Theorem [2l the Euler-Lagrange vector 
field is the unique semispray of type 1, Xc G X{T^Q), such that Xc o pr^ oj^ = Tpr^^ oX o j^. 
Thus this vector field Xc is locally given by 

d d 8 1,2 \ 9 

Xc = qiT, V 92^ V 93^ ('^ 90 + 92) ^ • 

dgo oqi ag2 7 c'93 

For the integral curves of Xc we know from Proposition [6] that if o" : R — )• W is an integral curve 
of X, then ac = P^i is an integral curve of Xc- Thus, if a is given locally by (j40p . then ac 
has the following local expression 

<ydt) = (9o(t),9i(i),92(t),93(t)) , (47) 



90 (i) 


= qi{t); 


91 (i) 


= 92 (t); 


92 (t) 


= 93(t); 


93 (t) 




p\t) 


= —uj'^q 


p\t) 


= 9i(i) 



and its components satisfy equations (l4Tj) . (1421) . (1431) and (l44]l . Notice that equations (1411) . (1421) 
and (|43p state that ac is the canonical lifting of a curve in the basis, that is, there exists a curve 
7: M — )• Q such that 7^ = ac- Furthermore, equation (|44p is the Euler-Lagrange equation for 
this system. 

Now, for the Hamiltonian solutions, as £ is a regular Lagrangian, Theorem [3] states that there 
exists a unique vector field Xh = TL^Xc G X(T*(TQ)) which is a solution to the Hamilton 
equation. Hence, it is given locally by 

9 9^91 d 

Xh = qi^ h 92^ qoTT-n + [qi-p ) TT- 

oqo oqi op^ ' op\ 

For the integral curves of X^-, Proposition [7] states that \i oc'- M — ?• T^Q is an integral curve of 
Xc coming from an integral curve cr of X, then fi/j = T L o ac is an integral curve of the vector 
field Xh- Therefore, if a is given locally by (j40]) . then ac is given by (|47|) and so a^ can be 
locally written 

cyhit) = {qo{t),qi{t),p''{t),p\t)) , 

and its components must satisfy equations (j4ip . (j42p . ()45p and (j46p . Notice that these equations 
are the standard Hamilton equations for this system. 
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4.2 The second-order relativistic particle 



Let us consider a relativistic particle whose action is proportional to its extrinsic curvature. This 
system was analyzed in [ 34^ [33| [H [31] , and here we study it using the Lagrangian-Hamiltonian 
unified formalism. 



The configuration space is a n-dimensional smooth manifold Q with local coordinates {q{^ 







1 ^ A ^ n. Then, if we take the natural set of coordinates on the higher-order tangent bundles 
over Q, the second-order Lagrangian function for this system, £ G C°°(T^Q), can be written 
locally as 



a 



a 



(48) 



{q[? '''''' ''''' ' {q\? 

where a is some nonzero constant. It is a singular Lagrangian, as we can see by computing the 
Hessian matrix of C with respect to g^, which is 



a 



[{{q\q^f - 2{q\f{qm q?qt + iq\fiq\q2)iq^qt 



2(gi)W 



qfq^) 



a 



-iq\?iq2fqiq^] 



[g - [qlfqtqf + mq\)qtq2 - iq\fq2q2' 



■\iB = A 



then after a long calculation we obtain that det 
regular Lagrangian. 



dq^dqi 



0. In particular, C is an almost- 



As this is a second-order dynamical system, the phase space that we consider is 



W = T3Qxtq T*(TQ) 




T*(TQ) 



As C is almost-regular, the "natural" phase space for this system would be T^Q Xtq -Poj where 
Po ^ T*{TQ) denotes the image of the Legendre-Ostrogradsky map. However, as we have a set 
of natural coordinates defined in W, it is easier to work in W and then to obtain the constraints 
as a consequence of the formalism. 

If uji G 17^(T*(T(5)) is the canonical symplectic form, we define the presymplectic form 
= pr2a;i G f2^(W), whose local expression is 

n = dqi)A dp^ + dq{ A dp} . 

The Hamiltonian function H £ C°°(W) is H = C — (p2°pi'i)*>C, where C is the coupling function, 
whose local expression is C (ggj Qij 021 q^iPi^Pi) ~ P^ql ~^Piqh then the Hamiltonian function 
can be written locally 



H (qlq\,qiqlplp}] 



pU\ + p}qi 



a 



{q\y 



[iqlf{ql?-iq\ql 



The dynamics for this system are described as the integral curves of vector fields X £ X(W) 
which are solutions to equation ([TH]). If we take a generic vector field X G X(W), given locally 
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by 



d 

9qo 



d 



+ F4 



d 



d 



d 



,0 



d 



taking into account that 
dH = qtdp^A + q2dp\ + 
+ 



a 



[{{<i,f{qh? - m^2?) qf + {Mmfq^' 



a 

-,i\2 



(qirVa 



{{q\fq^ - {qlq^)q^) 



dqi 



dq^ ; 



from the dynamical equation we obtain the following linear systems for the coefficients of X 

fo = qt (49) 



qt 

qt 




a 



a 



{q\yV9 



{{ql?q^ - {<fA)qt) 



- mqi?) qt + m2M?ql 
. 



(50) 
(51) 
(52) 

(53) 



Note that from equations (09]) and (|50p we obtain the condition of semispray of type 2 for X. 
Furthermore, equations (I53p are algebraic relations between the coordinates in W stating that 
the vector field X is defined along a submanifold Wo that is identified with the graph of the 
Legendre-Ostrogradsky map, as we stated in Propositions [2] and [3l Thus, the vector field X is 
given locally by 

ri . fi .ft .ft ^ fi 

(54) 



^ — qi TTT + q^ + + -f^s '^^r^ + lt^ 



dq^ ' dq^ ' dq^ ' dq^ ' ^^dp\ ' 

where the functions G\ are determined by (j52p . As we want to recover the Lagrangian solutions 
from the vector field X, we must require X to be a semispray of type 1. This condition reduces 
the set of vector fields X E X(>V) given by ([51]) to the following ones 



— qi tta + ^2 tta + ^3 + -^r-^ + ^ / 



dqi 



dqt 



dqi 



dqi 



dp\ 



Notice that the functions are not determinated until the tangency of the vector field X on Wo 
is required. Now, the Legendre-Ostrogradsky transformation is the map J-C: T^Q 
locally given by 



(55) 

iWo 

T*(TQ) 



FC*{p 



A) 



dC 



di 



dC 



dC 

dqi 



dr {p\) 



a 



FC*{p] 



(qiyva 

a 

(gi)VF 

dC a 



[{{q^?9 + {qi?iq^?{q\ql) - iq'ifi^iqDiqlql)) qi] + 



dqi {q\?^ 



mq\??{q\q\) - {q'lfiq'iqWiql) - {q\q\)9) qi - {q\?9qi] 
Wi?qi-{<iiq\)qi] , 



and, in fact, C is an almost-regular Lagrangian. Thus, from the expression in local coordinates of 
the map FC, we obtain the (primary) constraints that define the closed submanifold Pq = ImFjC, 
which are 



p}q\ = 



6f 



1\2 



iPt) 



a 



iq\) 



i \2 



; 



(56) 



P.D. Prieto-Martinez, N. Roman-Roy: Unified formalism for higher-order systems. 



31 



Pq. Then, the submanifold Wo = graph J^£o is defined by 
{p e W : CoiP) = ^Hp) = <l>f\p) = <P2\P) = 0, 1 ^ A ^ dimQ} 
Tli*-p\. The diagram for this situation is 




Notice that Wo is a submanifold of T'^Q xxq Pq^ and that Wo is the real phase space of the 
system, where the dynamics take place. 

Next we compute the tangency condition for X G j£(W) given locally by (fSSl) on the sub- 
manifold Wo ^ Wp„ ^ W, by checking if the following identities hold 

LWCo''|w„=0 , L(X)Ci^|w. = (57) 
L(X)(/.fVo=0 , L(X)4°V. = 0. (58) 

As we have seen in Section [3.2. 11 equations (j57p give us the Lagrangian equations for the vector 
field X. However, equations (I58p do not hold since 

L(X)<) = L(X)(pigl) = -pOgi , L(X)4°) = L(X)((pi)2 - aV(gi)') = -2ftV , 
and hence we obtain two first-generation secondary constraints 

tf^=P°9l = , <^f^v\p\ = ^ (59) 

that define a new submanifold Wi ^ Wq. Now, checking the tangency of the vector field X to 
this new submanifold, we obtain 

L(X)0i^)=L(X)(p,V) = O , L(X)<^«=L(X)(p0pi) = -(p0)2, 
and a second-generation secondary constraint appears 

0(2) ^ (p0)2 ^ ^ (go) 

which defines a new submanifold W2 ^ Wi. Finally, the tangency of the vector field X on this 
submanifold gives no new constraints, since 

L(X)0(2)=L(X)((pO)2)=O. 

So we have two primary constraints ()56p . two first-generation secondary constraints (j59p . and 
a single second-generation secondary constraint (I60p . Notice that these five constraints only 



Let TLo : T^Q ^ 

Wo = 

where = — 
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depend on qf-, p^^ and p\, and so they are pr2-projectable. Thus, we have the following diagram 



where 




Pi = [pePo: (^?{p) = ^?{p) = 0} = pr2(Wi 
P2 = {pePo: <t)^^\p) = 0} = pr2(W2) 

51 = Fc-\Pi) = wi{m) 

52 = TC-\P2) = wAy^2) . 



Focusing only on the Legendre-Ostrogradsky map, and ignoring the unified part of the diagram, 
we have 



rc 



T*(TQ) 




Notice that we still have to check (j57p . As we have seen in Section [3.2.11 we will obtain the 
following equations 



dq^dq^ dq^ 



dqi 



^1 '^^^d^ 



y {Fi^ - qi) dT 



dq^dq^ 



[Fi - qi) 



dq^dq^ 



= 
(61) 

(62) 
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As we have already required the vector field X to be a semispray of type 1, equations (|62|) are 
satisfied identically and equations ([6T]) become 



A long calculation shows that this equation is compatible and so no new constraints arise. Thus, 
we have no Lagrangian constraint appearing from the semispray condition. If some constraint 
had appeared, it would not be .FjCo-projectable (see 



Thus, the vector fields X G X(W) given locally by ([55]) which are solutions to the equation 

[i{X)n-dH\\y,^=Q , 

are tangent to the submanifold W2 ^ Wo- Therefore, taking the vector fields Xq € X(W2) such 
that Tj2 o Xq = X o j2, the form = (jWp^ ° jo ° ji ° J2)*^j and the canonical Hamiltonian 
function Ho = {jwp^ ° jo ° ji ° j2)*H, the above equation leads to 

i{Xo)no -dHo = o, 

but a simple calculation in local coordinates shows that Ho = 0, and thus the last equation 
becomes i{Xo)^lo = 0. 

One can easily check that, if the semispray condition is not required at the beginning and 
we perform all this procedure with the vector field given by (|54p . the final result is the same. 
This means that, in this case, the semispray condition does not give any additional constraint. 

As final results, we recover the Lagrangian and Hamiltonian vector fields from the vector 
field X G X{W). For the Lagrangian vector field, by using Lemma [3] and Theorem [2] we obtain 
a semispray of type 2, Xc G ^(T^Q), tangent to S2. Thus, requiring the condition of semispray 
of type 1 to be satisfied (perhaps on another submanifold M2 ^ ^2), the local expression for 
the vector field X^ is 

A d A d A d d 



dqo dqf dq^ dq 



3 



where the functions are determined by ()63p . For the Hamiltonian vector fields, recall that C 
is an almost-regular Lagrangian function. Thus, we know that there are Euler-Lagrange vector 
fields which are J-"£o-projectable on tangent to P2 and solutions to the Hamilton equation. 



5 Conclusions and outlook 

After introducing the natural geometric structures needed for describing higher-order autonomous 
dynamical systems, we review their Lagrangian and Hamiltonian formalisms, following the ex- 
position made in [l7j. 

The main contribution of this work is that we develop the Lagrangian-Hamiltonian unified 
formalism for higher-order dynamical systems, following the ideas of the original article [39]. We 
pay special attention to showing how the Lagrangian and Hamiltonian dynamics are recovered 
from this, both for regular and singular systems. 

A first consideration is to discuss the fundamental differences between the first-order and 
the higher-order unified Lagrangian-Hamiltonian formalisms. In particular: 
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• As there is no canonical pairing between the elements of T^'^^^Q of T*(T'^^^Q), in 
order to define the higher-order coupling function C in an intrinsic way, we use the canonical 
injection that transforms a point in T^^^^Q into a tangent vector along T''~^Q. 

• When the equations that define the Legendre-Ostrogradsky map are recovered from the 
unified formalism (both in the characterization of the compatibility submanifold Wo as 
the graph of J^C, and in the equations in local coordinates of the vector field X G X{W) 
solution to the dynamical equations), the only equations that are recovered are those that 
define the highest order momentum coordinates, and the remaining equations that define 
the map must be recovered using the relations between the momentum coordinates. 

• The regularity of the Lagrangian function is more relevant in the higher-order case, be- 
cause the condition of semispray of type 1 (the holonomy condition) of the Lagrangian 
vector field cannot be deduced from the dynamical equations if the Lagrangian is singular, 
unlike the first-order case, where this holonomy condition is deduced straightforwardly 
from the equations independently of the regularity of the Lagrangian function. When the 
Lagrangian is singular, we can only ensure that the Lagrangian vector field is a semispray 
of type k. It is therefore necessary, in general, to require the condition of semispray of 
type 1 as an additional condition. 

Then, for regular Lagrangian systems, when the tangency condition of the vector field 
X G X(yV) solution in the unified formalism along the submanifold Wo is required, we 
obtain not only the Euler-Lagrange equations for the vector field, but also the remaining 
k — 1 systems of equations that the vector field must satisfy to be a semispray of type 1. 

As we point out in the introduction, a previous and quick presentation of a unified formalism 
for higher-order systems was outlined in [13]. Our formalism differs from this one, since in that 
article the authors take T^Q (B^k-ig T*{T''~^Q) as the phase space in the unified formalism, 
instead of ours, which is T^'^-iQ T*(T'^ ^Q). This is a significant difference, since 

when we want to recover the dynamical solutions of the Lagrangian formalism from the unified 
formalism, the Lagrangian phase space is T'^^^^Q, instead of T^Q, which is the bundle where 
the Lagrangian function is defined. This fact makes it more natural to obtain the Lagrangian 
dynamics as well as the Hamiltonian dynamics, which in turn is obtained from the Lagrangian 
one using the Legendre-Ostrogradsky map. 

By using any suitable generalization of some of the several formalisms for first-order non- 
autonomous dynamical systems [U El [19] , a future avenue of research consists in generalizing 
this unified formalism for higher-order non-autonomous dynamical systems. This generalization 
should also be recovered as a particular case of the corresponding unified formalism for higher- 
order classical field theories. As regards this topic, a proposal for a unified formalism for higher- 
order classical field theories has recently been made [101 I41j . which is based on the model 
presented in [13] . This formulation allows us to improve some previous models for describing the 
Lagrangian and Hamiltonian formalisms. Nevertheless, some ambiguities arise when considering 
the solutions of the field equations. We hope that a suitable extension of our formalism to field 
theories will enable these difficulties to be overcome and complete the model given in [101 [41] . 
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